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Abstract. The deuteron-proton elastic scattering is studied in the multiple scattering expansion formalism. 
The contributions of the one-nucleon-exchange, single- and double scattering are taken into account. The 
Love and Franey parameterization of the nucleon-nucleon t-matrix is used, that gives an opportunity 
to include the off-energy-shell effects into calculations. Differential cross sections are considered at four 
energies, Td = 390, 500, 880, 1200 MeV. The obtained results are compared with the experimental data. 
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1 Introduction 

The study of the deuteron-proton elastic scattering has 
a longtime story. The first nucleon-deuteron experiments 
were performed already in fifties of the previous century 
[T]-[7]. Differential cross sections and polarization 

[S]-[7] were measured at few hundred MeV. Nowadays this 
reaction is still the subject of investigations [8]- [10]. This 
process is the simplest example of the hadron nucleus col- 
lision that is why the interest to this reaction is justified. A 
number of experiments on deuteron- nucleon elastic scat- 
tering is aimed at getting some information about the 
deuteron wave function and nucleon-nucleon amplitudes 
from Nd scattering observables. Moreover, the study of 
the reaction mechanisms, investigations of the few-body 
scattering dynamics are also very important to understand 
the nature of nuclear interactions. 

A good theoretical description of the deuteron- nucleon 
process was obtained for low energies, where the multiple 
scattering formalism based on the solution of the Faddeev 
equations, has been applied to solve this problem [TT] . 
However, at the energies above 150 MeV there is some dis- 
crepancy between the experimental data and theoretical 
predictions in the minimum of the differential cross sec- 
tion. At present many efforts are undertaken to extend the 
Faddeev calculation technique into the relativistic regime 
[12] -[H]. But up to now there are no reasonable descrip- 
tions of the experimental data obtained in the Faddeev 
equation framework at intermediate energies. 

Also p-d scattering at low energies was considered in 
the approach based on the solution of the three-particle 
Schrodinger equation using the Kohn variational principle 
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(KVP)[l6].[17j. Special attention in these works was given 
to the study of the Coulomb effects. It has been shown 
that at the energies below 30 MeV the influence of the 
Coulomb interaction is appreciable, while it considerably 
reduces at T^b = 65 MeV [T7] . 

The high-energy deuteron-proton scattering in the for- 
ward hemisphere is successfully described by the Glauber 
theory which takes both single and double interactions 
into account [TH] , [TH] • In [5D] it is shown that filling of the 
minimum, due to the interference between the single- and 
double-scattering amplitudes, is explained by the presence 
of the D-wave in the deuteron wave function. 

However, as it is shown in [21], [22], the off-energy- 
shell effects begin to play an important role at scattering 
angles larger than 30" in cm. In PSI the deuteron-proton 
backward elastic scattering was considered in the Bethe- 
Salpeter approach. Here various relativistic effects were 
studied. It has been shown, that the relativistic correc- 
tions coming from the negative energy P-waves are negli- 
gible while Lorentz boost effects become evident already 
at Piab ^ 0.2-^0.3 GeV/c (here,F;afc is the scattered proton 
momentum in the laboratory frame). 

The present paper considers the intermediate energy 
range from 200 MeV up to 600 MeV of the initial nucleon. 
On the one hand, the energies are not large enough to ap- 
ply the Glauber theory. On the other hand, in this region it 
is already necessary to take relativistic effects into account 
and as a consequence the Faddeev calculation technique 
does not work properly at these energies. Nevertheless, at 
these energies it is still possible to use a non-relativistic 
deuteron wave function due to transfer into the deuteron 
Breit frame. Here, it is also very important to describe 
correctly the nucleon- nucleon vertices, namely the spin 
structure as well as angular and energy dependences. 
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In the previous paper [22j the polarization observables, 
such as vector and tensor analyzing powers, were studied 
in the presented approach. In this work only the differ- 
ential cross sections are considered at four energies. The 
energy range has been chosen to demonstrate the region, 
where this model is justified. 

Sect. 2 gives the expression for the differential cross sec- 
tion and general kinematical definitions. The theoretical 
formalism is presented in Sect. 3. We consider a multiply- 
scattering series up to the second-order terms of the nucleon- 
nucleon f-matrix. The iterations of the Alt-Grassberger- 
Sandhas equations is performed to get this expansion. The 
wave function of the moving deuteron, which depends on 
the two variables, is used in the calculations. 

We apply the NN t-matrix constructed by Love and 
Franey [32! in order to describe the nucleon-nucleon in- 
teractions. This parameterization allows one to take the 
energy and angular dependences into account. Moreover, 
it is possible to extend this approach to the off-energy- 
shell region. The results of the calculations are discussed 
in Sect. 4. The conclusions are given in Sect. 5. 



2 Kinematics 

A cross section of the dp-scattering is expressed via the 
squared reaction amplitude: 



3 General formalism 

According to the three-body collision theory, the ampli- 
tude of the deuteron-proton elastic scattering J is defined 
by the matrix element of the transition operator Un: 

Udp^dp = 5{Ed + Ep — E'^ — E'p)J = 

< 1(23)|[1 - Pi2 - Pi3]f/ii|l(23) > . (6) 

Here, the state |1(23) > corresponds to the configuration, 
when nucleons 2 and 3 form the deuteron state and nu- 
cleon 1 is free. Emergence of the permutation operators 
for two nucleons Pij reflects the fact that the initial and 
final states are antisymmetric due to the two particles ex- 
change. 

The transition operators for rearrangement scatter- 
ing are defined by the Alt-Grassberger-Sandhas equations 



s{Pd+p-p'd-p'y 



E'^E'p J EdEp 



E'p) 



(1) 



This expression is valid for an arbitrary frame. Here {Ed, Pd), 
{Ep,p) are the energies and momenta of the initial deuteron 
and proton, respectively. The corresponding variables for 
the final particles are defined with the primed symbols. 
Starting from Eq. ([T|) one can obtain a formula for the dif- 
ferential cross section in the center-of-mass, 



da 



(27r) 



1 1 



E'^E^ J y/ EdEp\ 



dQ* " "' 6 s 
It should be emphasized that expression 



(2) 



I E'^E'p J \J EdEp^ is invariant. Since the reaction am- 
plitude J is calculated in the Breit frame, the energies 
should be also defined in the same system: 



Ed 



E'd 



-t/4, 



Ep - E'p 



{pQ)= Q' 



(3) 
(4) 



where Aid and mjv correspond to the deuteron and nu- 
cleon masses, respectively, and Q, p are momenta of the 
deuteron and proton in the Breit frame. The Mandel- 
stam variables s and t are defined through the laboratory 
deuteron energy and center-of-mass angle, 9*, 



s = {Pd+pf = Ml 
t={Pd-P'df = -- 



m 



N ^ 
2 

N t T>lab\2 



2mNE'f 



(Prr(l-cosr). (5) 



Uii ^ t2gQU2i + hgoUsi, 

U21 = 9q^ + hgoUii + tsgoUsi, 



(7) 



where ti = t{2, 3), etc., is the t-matrix of the two-nucleon 
interaction and go is the free three-particle propagator. 
The indices ij for the transition operators Uij denote free 
particles i and j in the final and initial states, respectively. 

Iterating these equations up to the ti-second-order terms, 
we can present the reaction amplitude as a sum of the fol- 
lowing three contributions: one-nucleon exchange, single 
scattering, and double scattering: 



Jdp^dp — JoNE + JsS + JdS, 

JoNE = -2 < l(23)|Pi25(ri|l(23) >, 
:7ss = 2<l(23)|ir'|l(23)>, 
JDS = 2<l(23)|irWr'|l(23)>, 



(8) 



where we have introduced notations for antisymmetrized 
operators tj^™ = [1— Pi3]t2 andfg^™ = [1— Pi2]i3- Schemat- 
ically, this sequence is shown in Fig.l. 

The convergence of the multiple scattering series was 
studied in ref . [14] , [15] . A reasonable agreement between 
the results of the calculations taking the second order cor- 
rections into account and the full Faddeev calculations, 
was obtained at the laboratory energies up to 2 GeV at 
least for the forward scattering angles. The analogous con- 
clusion can be done for the differential cross section of the 
elastic scattering at Eiab — 400 MeV for the scattering an- 
gle up to 140° [14] . It gives us the ground to suppose that 
the reduced series ([U is a good approximation to describe 
dp-elastic scattering at the energies above 400 MeV. 

After straightforward calculations we have got the fol- 
lowing expressions: 

for the one-nucleon-exchange amplitude - 



(9) 



JoNE - -\iEd -Ep- ^Jnij^+p^-Q"^) ■ 
< p'^i'; ~QM'd\n\{23)[l + {(7i(T2)]nd{23)\QMd;Pti > , 
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Fig. 1. The diagrams are taken into consideration: one-nucleon-exchange (a), single scattering (b), and double scattering (c) 
graphs. 



for the single scattering amplitude - 

Jss = Jdq' < -QM'a\f2l\q ^2, Q q '/i3 > (10) 
< p V, Q q '^i'3\^tl^{E) + ^t%iE)\pti,Q- q ^3 > 

<q'ti2,Q~q'i^3\i^d\QMd> . 

for the double scattering amplitude - 

Jds^ Jdq' J dq" < -QM'jn\ (11) 

I - Q - q V2,g " X p' i^':~Q - q " i^,q" M3I 

ti,{E')t\,{E) + [t\,{E') + tUE')][tUE) + t%{E)]/4 

Ed + Ep-Ei-E2-E'^ + le 
\p fj,,q' fi2,Q- g Vs > 
< q ' //2, Q q' tJi.'i\f^d\QMd > ■ 

Here the spin projections are denoted as fi for nucleons 
and A4(i ~ for deuterons. Henceforth, all summations over 



dummy discrete indices are implied. The superscript of the 
t-matrix corresponds to the isotopic spin of the nucleon- 
nucleon state. The argument of the A^iV-matrix is defined 
as the three-nucleon on-shell energy excluding the energy 
of the nucleon which does not participate in the interac- 
tion: 

E — Ed + Ep — E2, E' = Ed + Ep — E'^. (12) 

The three-nucleon free propagator in Eq. (|lip can be de- 
composed on two terms using the well-known formula: 



Ed + Ep- El- E2- E!^+ie ' ' 

V , 

Ed + Ep — El — E2 — i?3 

-iTrS{Ed + Ep-Ei-E2- E'^). 

The principal value part is often neglected to simplify 
the further calculations. However, as it was shown in [22], 
it is very important to take into account the full represen- 
tation of the three-nucleon propagator, especially, at the 
scattering angles above 40°. 
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All the calculations are performed in the deuteron Breit 
frame, where the deuterons move in opposite directions 
with equal momenta (Fig.l). It allows us to minimize the 
relative momenta of the nucleons in the both deuterons. 
As a consequence, the non-relativistic deuteron wave func- 
tion can be applied in the energy range under considera- 
tion. 

In the rest frame the non-relativistic wave function of 
the deuteron depends only on one variable po, which is the 
relative momentum of the outgoing proton and neutron: 

< flpfin\i^d\Md >= -T= < ^J'p^J■n\{u{po) + (14) 
\/47r 



wjpo 
V8 



■[3(cripo)(o'2Po) - (cTi cr2)]}\Md >, 



where u{po) and w(po) describe the 5" and D components 
of the deuteron wave function [26j , [27] , [28] , po is the unit 
vector in po direction. 

In order to get the wave function of the moving deuteron, 
it is necessary to apply the Lorenz transformations for the 
kinematical variables and Wigner rotations for the spin 
states. This procedure has been expounded in [22j, and, 
here, we give only the result: 

< pi m,P2^J■2\^^d\Q,Md >=< Pi Hi,P2H2\gi{k,Q) + 
g2{k, Q){eTin){(T2n) + g3{k, Q){(t^ct2.) + (15 
+ g4{k, Q){(Tik){(T2k) + g5{k, Q){{(Ti + (72)n) + 
+ geik,Q)[{(Tik){cT2,n x k) + {(Ti,n x k){cr2k)]\Q,Md > 

Note, the wave function of the moving deuteron is the 
function of two variables: the deuteron momentum Q and 
neutron-proton relative momentum k - 



The amplitudes A,B,C,D,F are the functions from the 
center-of-mass energy and scattering angle. The radial 
parts of these amplitudes are taken as a sum of Yukawa 
terms. A new fit of the model parameters was done [33] 
according to the current phase-shift-analysis data SP07 

Since the matrix elements are expressed through the ef- 
fective iVA^-interaction operators sandwiched between the 
initial and final plane-wave states, this construction can 
be extended to the off-shell case allowing the initial and 
final states to get the current values of the k and n' . Obvi- 
ously, this extrapolation does not change the general spin 
structure. 

To use definition (jl7[) for the nucleon-nucleon matrix in 
the further calculation, it is necessary to relate i-matrices 
in the two-nucleons cm. and in the deuteron Breit frame. 
This relation can be obtained by using the Lorenz trans- 
formation and Wigner rotation technique as it was done 
for the deuteron wave function. This problem has been 
considered in detail in |35] , [22| . 



< p'p'^]^l'^i'r^\t{E)\pp3]^l^lz >=Nt < K'^I[^i'2\ 

Wl^{p')Wl^M)ic.m.{V^^) Vt^l/2(p)W^l/2(P3)kMlM2 > ■ 

(19) 

ilere, momentum notations correspond to the ones given 
in Fig. lb. The Wigner rotations of the initial and final 
states are performed around n and n' axes, respectively: 



n 



\P X Pal 



n = 



P X P3 
|p ' X P3 'I 



(20) 



{E2 



k = 



P2, 

F V^/2)pi - {E, 



The Wigner rotation operator in the spin space of the z-th 
nucleon has the standard form: 



%/^/2)p2 



Ei+E2 + ^/^ 



(16) 



Here (i?i,pi) and (i?2,P2) are the energy and momentum 
of the neutron and proton in the deuteron, and s is the 
squared sum of the neutron and proton 4-momenta, s = 
[pi +P2)^. The normal to (pi,p2)-plane is denoted as n. 

In general, functions gi can be obtained by solving a 
relativistic equation, as it was done in |29j-[31j in the 
light front model. But in the present paper a usual non- 
relativistic DWF is taken as an input. Therefore, functions 
Qi are defined as the linear combinations of u and w {S- 
and Z^-waves). 

In order to describe the nucleon- nucleon interaction in 
a wide energy region, we have used the Love and Franey 
parameterization [32' of the A^A^ t-matrix defined in the 
center-of-mass as 

< K'fJ.[fl'2\tc.m.\l^^J■l^J■2 > = < k'^\^'^\A -I- (17) 
B{a^N*){(T2N*) + C(cri + era) • N* + D{(T:i_q*){a2q*) + 

i^((TiO*)(<T2Q*)|K/ilM2 > ■ 

Here the orthonormal basis is combinations of the nucle- 
ons' relative momenta in the initial k and final k' states: 



Wi/2{Pi,u) = exp{-iWi(niCri)/2} 
cos(wi/2)[l - ilniCTi)tg{uji/2)], 



(21) 



where uji is the angle of this rotation. The normalization 
factor M provides conditions of orthonormality and com- 
pleteness and is defined by the Jacobian of the transfor- 
mations. The kinematical factor T is due to the relation 
between two off-energy shell t-matrices, one of them de- 
pends on zero total momentum |36| , [22] . 

The argument of the t-matrix in the right-hand side of 
Eg. p^ is the square root of the Mandelstam variable s, 
which is defined through the kinematical variables in the 
deuteron Breit frame taken on the mass-shell: 



E' 



(22) 



K- k' _ K + k' 



^ /v X txi , ^ , 

N* = -. (18) 



where iiT = p -I- pa is the total momentum of the nucleon- 
nucleon pair in the frame of the calculation. 



4 Results and discussion 

The results of the calculations for the differential cross 
sections are presented in Figs(2][5j The four deuteron ki- 
netic energies, 390, 500, 880 and 1200 MeV, have been 
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0* 

Fig. 2. The differential cross section at the deuteron kinetic 
energy of 390 MeV as a function of the cm. scattering angle. 
The dashed line corresponds to the calculations including only 
one-nucleon-exchange and single scattering diagrams into con- 
sideration. The solid and dashed-dotted lines are the results of 
the calculations taking into account also the double scattering 
diagram with the principal value part of the free propagator 
and without it, respectively. The data are taken from [9] (■) 
and [37] (O)- 



considered. All calculations were performed with the CD- 
Bonn deuteron wave function j27j. The dashed line in 
FigsHHS] corresponds to the calculations taking into ac- 
count only ONE and single-scattering contributions. The 
dash-dotted line represents the calculation including the 
double-scattering term without the principal value part in 
the three-nucleon propagator, and the solid line stands for 
the calculations with the full three-nucleon propagator. 

One can see that the results taking the double scatter- 
ing into account and without it, are close to each other 
up to 50-60 degrees of the scattering angle for all the 
considered energies. At these angles the single-scattering 
gives the main contribution to the differential cross sec- 
tion. Then the difference between these results increases. 
It should be noted that the double scattering contribution 
is not large enough for the deuteron energy of 390 MeV 
(Fig. 2). As a consequence, the curves corresponding to the 
results with and without the principal value part of the 
three-nucleon propagator, are practically undistinguished. 
However, the double-scattering effect increases with the 




0* 

Fig. 3. The differential cross section at the deuteron kinetic 
energy of 500 MeV as a function of the cm. scattering angle. 
The lines are the same as in Fig[2] The data are taken from[8]. 

energy. Already for the energy of 500 MeV (Fig. 3) the dif- 
ference between the results taking DS into consideration 
and without it, is significant - it is about 2 times in min- 
imum of the differential cross section. For higher energies 
this factor is even larger: it is about 10 times for Td = 880 
MeV (Fig. 4), and about 20 times for the energy equal to 
1200 MeV (Fig.5). 

The difference between the results taking into account 
the principal value part of the three-nucleon propagator 
("full" calculation) and without it, also increases with the 
energy. At small angles this deviation is practically invis- 
ible, while it is significant for the scattering in the back- 
ward hemisphere, 9* > 90°, where the double scattering 
effect plays the main role. Such a discrepancy can be, 
probably, explained by the relativistic effects related with 
the off-energy-shell behaviour of the nucleon-nucleon t- 
matrices, and Fermi motion of the nucleons in the deuterons. 
It is interesting to emphasize that the curves correspond- 
ing to the calculations without the principal value part of 
the free propagator are close to the ONE-|-SS-results at 
the scattering angle above 160" at all the energies under 
consideration. 

The deviation of the theoretical predictions from the 
experimental data are observed at the scattering angles 
above 140°. One can assume some reasons of this distinc- 
tion. It is possible that the multiple scattering series is not 
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Fig. 4. The differential cross section at the deuteron kinetic 
energy of 880 MeV as a function of the cm. scattering angle. 
The lines are the same as in Fig[2] The data are taken from 
[38] (•), and[39](A). 



Fig. 5. The differential cross section at the deuteron kinetic 
energy of 1200 MeV as a function of the cm. scattering angle. 
The lines are the same as in Fig[2] The data are taken from [40]. 



yet converged in the backward direction, and it is neces- 
sary to take the more orders of this series into account. 
This is supposed according to the results of the investiga- 
tion performed in ref.[14j. 

The agreement between the experimental data and 
theoretical predictions could be also improved, if the three- 
nucleon forces are included into consideration. In fact, a 
good description of the differential cross sections, both 
in the diffraction minimum and the backward direction, 
was obtained at the energies below 400 MeV [41], when 
the three-nucleon forces were taken into account. Unfortu- 
nately, at present there are no such calculations at higher 
energies. However, it should be noted, that the three- 
nucleon forces can be effectively included through the A- 
isobar excitation. The employment of this approach [32] 
gives the results, which are in agreement with the re- 
sults obtained in [41]. The essential contribution of the 
Zi-isobar excitation into the differential cross section was 
also demonstrated in ref. [43] , where the pd-backward elas- 
tic scattering was studied at the high energies. 

Nevertheless, the results of the "full" calculations are 
in a reasonable agreement with the experimental data up 
to the scattering angle of 140° at all the considered ener- 
gies. 



5 Conclusion 



In the paper we have presented a method to calculate 
the amplitude of the deuteron-proton elastic scattering 
at intermediate energies. Special attention was given to 
the questions connected with the relativistic effects. The 
transformation of the deuteron wave function in the rest 
frame to a moving system was performed, that allowed us 
to use the non-relativistic DWF at rather high energies. In 
order to describe nucleon-nucleon interactions in a wide 
energy range, we have used parameterization of the NN 
t-matrix. The spin transformation technique has been also 
applied to relate this t-matrix given in the cm. to that in 
the reference frame. 

Using this method we have managed to describe the 
experimental data on the differential cross section at four 
energies: 395, 500, 880, and 1200 MeV. Good agreements 
have been obtained between the experimental data and 
the theoretical model calculations taking into considera- 
tion the one-nucleon-exchange, single- and double-scattering 
for all the four energies up to the scattering angle of 140°. 
The distinctions between the data and theoretical pre- 
dictions at the backward angles should be studied more 
precisely, that is the subject for further investigations. 
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